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-—L ' Abstract. We present a model for the OH megamaser emission in the starburst 

^— «, , galaxy IIIZw35. The observed diffuse and compact OH maser components in this 

^si ■ source are explained by a single phase of unsaturated clumpy gas distributed in 

a thin ring structure and amplifying background continuum. We emphasize the 
importance of dumpiness in the OH masing medium, an effect that has not been 
fully appreciated previously. 

The model explains why multiple bright spots are seen only at the ring tangents 

while smoother emission is found elsewhere. Both the observed velocity gradients 

and the line to continuum ratios around the ring enquire a geometry where most of 

the seed photons come from a continuum emission which lies outside the OH ring. To 

►^ . explain both the OH and continuum brightness, free-free absorbing gas is required 

t***^ ' along the ring axis to partially absorb the far side of the ring. It is proposed that the 

l/^ required geometry arises from an inwardly propagating ring of starburst activity. 

^' 

\^ • Keywords: starburst, maser 

o 
o 

T^ • 1. Introduction 

O I Extragalactic OH megamasers usually present both extended diffuse 

-4—' ' emission and very compact and strong features. One of the clearest 

C^ , cases of an OH megamaser showing both compact and diffuse masers 

is niZw35. Two groups of compact masers were detected in VLBI 
(Trotter et al. 1998; Diamond et al. 1989) recovering nearly half of the 
r> I total emission seen at MERLIN scales (Montgomery and Cohen 1992). 

C^ ■ Aiming to determine the location of the missing diffuse component 

MERLIN+EVN observations were conducted on this source (Pihlstrom 
et al. 2001). These observations revealed the best example of a rotating 
OH maser ring that has yet been found. 

It has been claimed that the two-phase characteristic of the maser 
features observed toward IIIZw35 is due to two coexisting phases of 
the OH gas (Diamond et al. 1989), one producing the smooth emission 
and the other responsible for the very strong and compact features. 
According to this hypothesis, and making the usual assumption that 
the diffuse component maser is operating in the unsaturated regime, 
the Line to Continuum Ratio (LCR) of each phase is exponentially 
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EVN+MERLIN 

1667 MHz OH inaset emission 




Figure 1. Left: Velocity averaged OH emission as seen at EVN+MERLIN resolution 
(Pihstrom et al. 2001). The contours are -1, 1, 2, 4 and 8 times the 3(7 rms noise of 2.7 
mjy/beam. Right: The proposed geometry. The OH gas is concentrated within the 
dark gray ring and most of the continuum is emitted within the surrounding light 
gray ring. The medium gray cones indicate the region where free-free absorption 
occurs. The units are normalized to the outer OH radius Ro. 



related to the path length through the amplifying gas. In the particular 
case of the smooth component in IIIZw35, the LCR in the east/west 
and north/south regions has been observed to be around 10 and 1000 
respectively. Additionally, the ratio of the path lengths between these 
regions is estimated to be approximately 8 therefore, the ratio of their 
LCRs should be in the order of 10^. Instead in order to match the 
observations, a clumpy medium is proposed so the column density of 
amplifying gas is proportional to the number of overlapped clumps 
along the line of sight. This number is the result of a random process 
that yields a mean which is linearly related to the projected path length 
through the gas. 



2. Description of the model 



2.1. Geometry 



We propose that the conditions for the OH gas to produce maser ampli- 
fication are satisfied within a thin ring structure with inner and outer 
radii Rj and Rq. Surrounding the OH, there is another concentric thin 
ring of radius Re where most of the continuum emission is produced. 
In addition, free-free absorption cones along the axis are assumed to 
account for the east-west asymmetry in the brightness distribution. 
This situation is depicted in Figure 1 and the observed geometrical 
parameters are summarized in Table I. 
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Table I. Observed geometrical parameters. 



Parameter 


Value 


Reference 


AR 


3Pc 


Diamond et al. 1998; Trotter et al. 1989 


AH 


6 Pc 


Diamond et al. 1998; Trotter et al. 1989 


Ri 


22 Pc 


Trotter et al. 1989 


Ro 


25 Pc 


Pihlstrom et al. 2001 


Re 


50 Pc 


Pihlstrom et al. 2001 


e 


60 Deg 


Pihlstrom et al. 2001 



The fact that the northern and southern compact masers seem to 
be elongated toward the east (Diamond et ah 1998; Trotter et ah 1989) 
is explained in our model by the larger amount of continuum behind 
the OH in these regions. Accordingly, the apparent east-west extention 
of the feature corresponds to half of the ring width AH whereas the 
north-south extention matches the difference between the outer and 
inner radii AR. The angle 6 between the axis of symmetry of the disk 
and the direction of the observer is estimated from the aspect ratio of 
the projected ring (Pihlstrom et al. 2001). 

2.2. Clumpy medium analysis 

The dumpiness of the medium is defined by the parameter n which 
defines the average number of clumps per cloud velocity width found 
in the radial direction within the ring. The average number of clumps 
in the direction of the observer is given by A(x,y) = n ■ L{x,y)/AR 
where L{x,y) is the projected path length through the OH gas. For 
any given line of sight, the number of observed spatial over lappings per 
velocity range is drawn from a Poisson random process of parameter A 
with Probability Mass Function Px (n) . If all the clumps have the same 
optical depth tq, then the expected amplification G along the line is 
given by the sum: 

N 

G=^PA(n)-e-° (1) 

n=0 

which in the limit of iV — > oo converges to e^^'^^''~^' (See Conway et 
al. in this volume). This result is inappropriate in our context since 
the maximum number of aligned clumps N is upper bounded either by 
geometrical or statistical constrains. For these reasons, the mean does 
not constitute a good estimator for the observations. A better choice 
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could be to use the median instead. A simplified approximation is to 
define an expected maximum N in order to truncate the sum in (1). 
Given a set of M independent realizations of the process, corresponding 
to M independent lines of sight within the beam, and requiring at least 
one occurrence of A^ overlaps to be observed, the expression for the 
maximum N becomes: 



iVA(M) = P-(l--i) 



Where P^ (x) is the Inverse Poisson Cumulative Distribution of pa- 
rameter A, and M is given by the ratio of the areas of the EVN+MERLIN 
to the VLBI synthesized beams when the clumps are assumed to be just 
resolved at the latter resolution. For the proposed geometry we define 
the following quantities and their observed values: 

Gi = J2 ^Ai (n) ■ e"""" « 10 Gain in Region 1 

n=0 

N2 

G2 = Yl ^A2(n) • e"""" w 50 Gain in Region 2 (3) 

n=0 

Gp = e^2To > 5Q0 Peak gain in Region 2 

Where Ni and N2 are the maximum number of cloud overlaps in regions 
1 and 2 given by (2) and A2 = i^ • Ai with K being the ratio of the 
averaged path lengths within both regions. Each of these equations 
define a locus in the Ai — tq plane where the observed values are 
matched. The intersection of those loci corresponds to the possible 
combinations of Ai and tq satisfying all three conditions simultaneously. 
Solving numerically, the solution is found to be Ai = 0.67 and tq = 1.69 
which corresponds with a mean free path between clumps i ~ 1.72 pc. 



3. Simulations 

Using the parameters obtained in the previous section as input for 
Monte Carlo simulations, artificial spectra and brightness distributions 
were produced. Random generators were set up to reproduce the spatial 
distribution of the individual clumps and the pseudo-turbulent velocity 
field. The resulting XYZ cube was transformed into a XYV cube using 
the cascade model shown in Figure 3. The results are summarized in 
Figures 4 and 5. 
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Figure 2. Left: Diagram of the model. The boxes indicate the positions of regions 1 
and 2. The outer radius of the continuum ring is not visible in this drawing. Right: 
Solutions of equations (3) shown in the Ai — tq space. The intersection of the dark 
and light gray lines is the solution of the system. 



i(x.y.v) 





Figure 3. Cascaded amplifiers model. The gain d of each amplifier is given by 
exp{ro-exp((«— «j)^/2(T^)} Where u = 10 km s~^ is the internal velocity dispersion of 
the clumps and Vj is the projected velocity of each clump constituted by a turbulent 
and a keplerian component, d are the continuum contributions between the clumps 
and Aff is the free-free absorption attenuation. All the parameters are functions of 



4. Conclusions 

The proposed geometry has managed to reproduce the observed bright- 
ness distribution and spectrum of the OH maser emission. The analyti- 
cal approach in section 2 has been useful in fine tuning the Monte Carlo 
simulations. However, some aspects of the observed velocity field need 
to be studied. Work is still in progress connecting the obtained results 
with the physics of masers and the developing of a consistent formalism 
for the analysis of clumpy media. 
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Figure 4- Integrated OH intensity I{x, y). The XYV cube was averaged along all the 
spectral channels. The boxes indicate the northern and southern regions which are 
shown enlarged in the right panel. Both regions present strong compact emission 
dominated by a few very bright spots. In addition, the arc shapes observed with 
VLBI are fairly well reproduced. The grayscale has been normalized to the brightest 
point in the south. 
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Figure 5. Left: Integrated OH intensity convolved down to EVN+MERLIN res- 
olution. The bridges of smooth emission connecting the north and the south are 
clearly visible as well as the east-west asymmetry produced by the free- free absorbing 
cones. The grayscale has been normalized to the brightest point in the south. Boxes 
indicate the positions where spectra were taken. Right: The spatial location and 
relative widths of the spectral features correspond to those observed by Pihlstrom 
et al. 2001. The vertical scale has been normalized to the peak in the south. 
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